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ABSTRACT 

Results of a study of the theoretically predicted and observed X-ray properties of 
local massive star clusters are presented, with a focus on understanding the mass and 
energy flow from these clusters into the ISM via a cluster wind. A simple theoretical 
model, based on the work of Chevalier & Clegg (1985), is used to predict the theoretical 
cluster properties, and these are compared to those obtained from recent Chandra 
observations. The model includes the effect of lower energy transfer efficiency and 
mass-loading. In spite of limited statistics, some general trends are indicated; the 
observed temperature of the diffuse X-ray emission is lower than that predicted from 
the stellar mass and energy input rates, but the predicted scaling of X-ray luminosity 
with cluster parameters is seen. The implications of these results are discussed. 

Key words: hydrodynamics: shock waves; open clusters and associations; stars: 
winds, outflows 



1 INTRODUCTION 

Super Star Clusters are dense clusters of young massive 
stars, first identified in NGC 1275 by Holtzman et al. (1992) 
using the Hubble Space Telescope (HST), and subsequently 
in a wide range of star-forming galaxies, such as merg- 
ing systems (NGC 4038/4039; Whitmore & Schweizer 1995), 
dwarf galaxies (Henize 2-10; Johnson et al. 2000), classical 
starbursts (M82; Gallagher & Smith 1999) amongst many 
other systems (see Whitmore 2000 for a review). These 
extragalactic star clusters can contain many thousands of 
very young, energetic stars, and have stellar densities far 
greater than those seen in normal OB associations. The 
values quoted by Whitmore (2000) indicate ages for these 
star clusters of typically 1 — lOMyr, radii typically in the 
range of ~ 1 — 6pc, total masses of stars in the cluster in 
the range 10 3 — 10 6 Mq ) , with the central stellar densities 
reaching up to ~ 10 5 Mq pc 3 . It is clear that in many galax- 
ies a substantial fraction of the ongoing star-formation (and 
hence mass/energy injection into the ISM) is occurring in 
SSCs (Origlia et al. 2001). The Galactic (or local analogues) 
of these extragalactic SSCs are objects such as NGC 3603, 
R136 in 30 Doradus and the Arches cluster near the Galactic 
Centre (see Figer, McLean & Morris 1999b). 

The component stars of an SSC are believed to be 
roughly coeval, and when the cluster is very young the mas- 
sive stars in the cluster will have strong stellar winds which 
will combine to produce an outflow from the cluster, which 
is referred to as the cluster wind. As the cluster ages mass 
and energy input from supernova explosions will begin to 



dominate (cf. Leitherer & Heckman 1995). In many star- 
burst galaxies, where a large number of SSCs are seen, this 
cluster wind will be an important contributor to energising 
the interstellar medium and possibly driving a hot super- 
wind or outflow, such as those seen in M82 or NGC 253 
(Lehnert, Heckman & Weaver 1999; Strickland et al. 2002). 
The efficiency with which stellar kinetic energy (from both 
stellar winds and supernovae) is converted into thermal en- 
ergy within the cluster, which in turn can drive an outflow, 
is an important and very difficult parameter to determine. 
In the literature, several different values or prescriptions for 
the thermalization efficiency have been used, often encapsu- 
lated in a parameter r/, which represents the fraction of the 
kinetic energy of stars and supernova in the cluster that is 
thermalized. The fraction of the kinetic energy that is not 
thermalized is assumed to be lost to the system, primarily 
via radiative losses. As we shall see radiative losses at X- 
ray energies are likely to very small, however, at UV and 
IR energies radiative losses are likely to be much larger. In 
fact, inhomogeneities in the complex flow within the clus- 
ter could easily lead to the formation of denser regions that 
could contribute to radiating energy away from the system. 

In the literature, Strickland & Stevens (2000) argued 
for very efficient thermalization (rj ~ 1; see also Chevalier & 
Clegg 1985), whereas Bradamante, Matteucci & D'Ercole 
(1998) argued for a much lower thermalization efficiency 
of a few per cent. Other, more complex, prescriptions of 
rj have also been proposed. For instance, Recchi, Matteucci 
& D'Ercole (2001) use a model with r\ being very low in 
the early evolution of a cluster, but as material is evacu- 
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ated by the cluster wind to create a low density region the 
thermalization efficiency rises to close to unity. In relation 
to LMC superbubbles, Oey (1996) noted that the bubble 
dynamics implied a relatively low thermalization efficiency 
(or put another way that there was less power driving the 
bubbles than implied by the stellar populations within the 
bubbles, see also Chu & Mac Low for a discussion of the 
X-ray properties of superbubbles). It would clearly be very 
useful to have more direct observational constraints on this 
process. 

Understanding the mass and energy loss processes from 
distant SSCs via the hot cluster wind can best be done at 
X-ray energies. However, because of their distance and the 
high stellar density, even Chandra has insufficient resolution 
to see what is going on in and around these systems, and 
resolve the diffuse emission from point source emission from 
stars in the cluster. However, local lower-mass analogues of 
these SSCs in our Galaxy or Local Group can be used to 
study the mass and energy loss processes from stellar clus- 
ters. By studying the diffuse X-ray properties of star clusters 
and the stellar energy injection rate from stellar populations 
in the cluster it should be possible to constrain the efficiency 
with which the clusters can convert kinetic energy to ther- 
mal energy. Although not perfect analogues, the values de- 
termined could have far-reaching consequences, not just for 
the understanding of the processes within the star clusters, 
but also for the energy emission from larger structures such 
as starburst galaxies, where many SSCs are present. 

It is also worth making the point that while the mass 
and energy injection rate from massive stars via their stellar 
winds can be estimated to a reasonable degree of accuracy, 
the same is less true for supernova. For instance, it is worth 
noting that the type Iln supernova SN1988Z is believed to 
have radiated ~ 10 52 erg in its early evolution (Aretxaga 
et al. 1999), see also Chevalier & Fransson 2001), which is 
rather more than the usually assumed values for the total 
energy injection from a SN. Consequently, looking at very 
young clusters, where SN injection does not dominate, has 
considerably more promise. 

In this paper the theory of cluster winds is devel- 
oped, particularly as it relates to X-ray emission from 
clusters (§ 2), and it is then applied to results from Chan- 
dra observations of nearby SSCs or smaller stellar clusters 
(§ 3). The results are discussed in § 4 and summarised in § 5. 



2 THE X-RAY EMISSION FROM A CLUSTER 
WIND 

When massive stars with strong winds are in close proxim- 
ity, such as in a stellar cluster, the winds collide both with 
each other and with the surrounding ISM, filling the sur- 
rounding volume with hot shocked gas. Eventually all of the 
available volume is filled and the hot gas escapes beyond 
the boundaries of the cluster. The development of such a 
cluster wind has been numerically modelled by Canto, Raga 
& Rodriguez (2000) in a simulation carried out assuming a 
spherical cluster consisting of ~ 30 stars (see also Ozernoy, 
Genzel & Usov 1997; Raga et al. 2001). 

The wind from a stellar cluster can also be described 
using the Chevalier & Clegg (1985) model originally used to 



describe outflows from starburst galaxies. The contrast be- 
tween these models is that Canto et al. (2000) also used hy- 
drodynamic calculations with mass and energy input from 
discrete stars within the cluster, whereas the Chevalier & 
Clegg (1985) model assumed mass and energy injection uni- 
formly throughout the starburst region. The dense star clus- 
ters considered here do not have a uniform distribution of 
stars and tend to be centrally concentrated, and mass and 
energy injection will tend to follow the stellar distribution. 
However, Canto et al. (2000), who did include a non-uniform 
stellar distribution showed that their numerical simulations 
generally reproduce the results of the Chevalier & Clegg type 
model, and the formulation of Chevalier & Clegg (1985) is 
adopted here. One important difference between the Canto 
et al. (2000) and Chevalier & Clegg (1985) models is that 
the Canto et al. (2000) have a much greater range of tem- 
peratures (and densities inside the cluster core). 

If mass and energy from stars is being injected via stel- 
lar winds into the volume of a stellar cluster of core radius 
R c then Chevalier & Clegg (1985) showed that the solution 
for the outflow from this region can be written as 
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for r > R c , where M is the Mach number of the flow, 7 is 
the adiabatic index (with 7 = 5/3 assumed hereafter), R c 
the radius of the star cluster and r is the radius from the 
center of star cluster. From these equations and the mass and 
energy injection rates, the velocity v(r), temperature T(r) 
and mass density p(r) can be determined. The technique 
to do so involves solving eqn.^or eqn.|5|as appropriate, for 
the Mach number at each radius, using the Newton-Raphson 
method. From the Mach number, and using the integrated 
forms of the mass and energy energy equations, the other 
variables can be determined (see Chevalier & Clegg 1985 for 
more details). Solar abundances and fully ionized material 
with a mean mass per particle of p, = 10~ 24 gm are also 
assumed. 

It is worth noting that the conditions in the central re- 
gions of the cluster are a good indicator of the recent history 
of the cluster. The flow time within the cluster is short, typ- 
ically a few 10 3 yr, whereas the characteristic timescale of 
the larger scale bubble being blown by the cluster is much 
longer (i.e. of order the age of the cluster). This means that 
the current X-ray properties should give a better handle on 
the transfer efficiency of stellar wind energy to the cluster 
wind. 



2.1 The Stellar Mass and Energy Injection Rate 

For the young star clusters under consideration here the 
energy injection is dominated by stellar winds, especially 
those of the massive OB and Wolf-Rayet (WR) stars, and 
we assume that no supernovae have as yet exploded and 
contributed to the cluster wind. At later epochs supernovae 
begin to dominate the energy injection (Leitherer & Heck- 
man 1995). 
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Figure 1. The radial structure for the cluster wind, as calculated from the Chevalier & Clegg (1985) model (see text for details). Shown 
are the radial profiles for the wind velocity (top left panel), gas temperature (top right), density (bottom left) and the cumulative X-ray 
luminosity (bottom right). The standard model (solid line) is for a star cluster with a stellar mass injection rate of M„ = 10~ 4 Mq yr _1 , 
V* = 2000 km s _1 , rj = 1 and no additional mass-loss (M co (d = 0). The second model is identical, except with r] = 0.1 (dashed line), 
and the third model (dotted line) has r\ = 1 and M ao id = 2 X 10~ 4 Mg yr _1 . 



In principal, from the current observed massive star 
population in the clusters, the mass and energy injection 
rate from the stellar winds of these stars can be estimated, 
and this is the procedure we adopt. However, as discussed 
elsewhere, there are some significant uncertainties in the de- 
termination of these values, related to uncertainties in un- 
derstanding the structure and evolution of massive stars. 

The bulk of the mass and energy injection is assumed 
to happen within a cluster core radius R c . If the cluster 
contains N* stars, each star with a mass-loss rate of M and 
wind terminal velocity Vi with i = 1 . . . N, then the total 
stellar wind mass injection rate, M», is 



M, = Mi 

i 

and the total kinetic energy injection rate is 



(3) 



(4) 



A mean weighted terminal velocity for the stars in the clus- 
ters is defined such that 



V, 



1/2 



(5) 



For some of the stellar clusters discussed in the next 
section the individual mass-loss rates and terminal velocities 
can be determined, and for those systems where there is 
insufficient observational evidence stellar evolution models, 
such as the Starburst99 code (Leitherer et al. 1999) can be 
used. 

Some qualifications should be made here. The method 
used for estimating the mass-loss rates and terminal veloc- 
ities of the cluster stars rests on our current best under- 
standing of the mass- loss parameters of massive stars. There 
remain significant deficiencies in our knowledge of the struc- 
ture and evolution of massive stars and their winds. For in- 
stance, the Starburst99 code only assumes a population of 
single stars and binary evolution could make a difference 
to the mass-loss rates and terminal velocity. Indeed, it is 
known that the characteristics of the winds of massive stars 
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evolve with time, and this snapshot could be missing impor- 
tant early phases of evolution when the winds were much 
stronger. In spite of these qualifications we shall use current 
best estimates and proceed accordingly. 

In a cluster although the optically visible stars (and 
their winds) will likely dominate the energy injection, these 
stars will not necessarily be the only sources of mass. Pro- 
tostars (or other optically obscured objects) will also con- 
tribute mass as material is ablated off them by the cluster 
wind or through radiation. Examples of such objects may 
be the ProPlyDs which have been seen in NGC 3603 and 
elsewhere (for instance, Brandner et al. 2000; Miicke et al. 
2002). This cold mass being injected into the cluster wind 
can be accounted for by a term M co id, so that the total 
mass-loss rate into the cluster is 

Mm = M, + M cold (6) 

This process can be considered as a mass-loading process 
(see Hartquist et al. 1997 and references therein). 

As discussed in the introduction, to account for the ef- 
fects of radiative losses in the conversion of stellar wind en- 
ergy into the cluster wind, we can introduce a parameter 
that we term the "thermalization efficiency" parameter rj, 
such that the energy available to drive the cluster wind is 
E t h = vE* , so that 
n, 

E th = V J2 \™ iV i = (7) 

i 

and allow rj to vary between and 1. This parameter allows 
for radiative losses. The energy radiated at X-ray energies 
is unlikely to be a major loss mechanism, but UV radiation, 
particularly from interfaces between hot and cold gas, may 
be a bigger contributor, as could IR radiation. As discussed 
by Recchi et al. (2001) we might expect this parameter to be 
relatively small in the early stages of cluster evolution, but 
to increase at later times. For simplicity, in these simulations 
r\ is assumed constant. The result of allowing rj to be less 
than unity and additional "cold" mass-injection will be to 
reduce the terminal cluster wind velocity below V*. This is 
turn increases the density and often leads to an enhanced 
X-ray luminosity. 

The assumption of a constant thermalization efficiency 
is also highly simplistic. Not only is the value likely to vary 
as a function of time, as the cluster wind evolves, it is also 
likely to vary as a function of location within the cluster, 
with a different thermalization efficiency in the denser inner 
regions to that in the more diffuse outer regions. Given the 
lack of detailed information as to how rj might vary, we shall 
treat it as a constant. So, subject to this recognition that this 
is a highly simplified treatment of a very complex problem 
we shall proceed. 

2.2 The Cluster X-ray Luminosity 

Observationally, at X-ray wavelengths, a young stellar clus- 
ter will have both point source emission associated with indi- 
vidual stars, such as from colliding winds in massive binaries, 
intrinsic emission from single early-type stars, emission asso- 
ciated with SNRs, X-ray binaries or pre-main sequence stars, 
and diffuse emission (due to the cluster wind). Here we con- 
centrate only on the diffuse X-ray emission, and only from 



younger clusters where we would not expect to see emission 
from SNRs and X-ray binaries. In the data presented later, 
the point sources have been removed from the X-ray data to 
leave only the diffuse emission. It is however possible that 
unresolved point sources, particularly from lower luminosity 
pre-main sequence (PMS) stars, will contaminate the diffuse 
emission. We will discuss this point more later. 

For a given model, the broad-band X-ray luminosity of 
the cluster will be given by 

L x = / 4:nr 2 n e riiA(T)dr (8) 
Jo 

where A(T) is the emissivity of gas at a temperature T, and 
n e and m are the electron and ion number density respec- 
tively. Given the solution for p(r), v(r) etc, eqn. [H]can be 
evaluated to derive the total X-ray luminosity (over a wide 
waveband) of the cluster. A more detailed analysis could 
calculate the luminosity in specific wavebands. However, the 
majority of the emission will fall in the Chandra waveband 
and we will not complicate matters further. 

In Fig.^the radial variation of the velocity and temper- 
ature for a stellar cluster are shown, where the total stellar 
mass-loss rate is 10 -4 M© yr _1 and V* = 2000 km s" 1 and 
the cluster core radius is lpc. These values are similar to 
those of the clusters discussed in this paper. In this model 
the central temperature is 5.8 x 10 7 K, the central ion num- 
ber density is 0.67cm~ 3 , and the X-ray luminosity of the 
region within a radius R c is Lx = 5.1 x 10 32 erg s -1 . Im- 
portantly for our purposes, the cluster itself dominates the 
X-ray emission and the volume outside R c contributes only 
15 per cent of the total X-ray luminosity. The reason for 
this is the sharp fall off in density (and temperature) outside 
the cluster, which more than compensates for the increasing 
volume of the emission region. For a simple analysis such as 
presented here this level of accuracy is sufficient. 

The central temperature of the cluster wind (To) in this 
model is given by 

=i.45xl0 7 ( ^ r V . (9) 

V K / ^ 1000 km s" 1 J V ; 

Note that the slightly different numerical constant as com- 
pared to Canto et al. (2000) is due to slightly different as- 
sumptions as to the mean mass per particle. 

From the radial solution the emission weighted tem- 
perature for the cluster region (T c ;) can also be calculated, 
which will be somewhat lower than the central temper- 
ature. For instance, for the standard model with V* = 
2000 km s _1 , the emission weighted temperature is T c i = 
5.1 x 10 7 K. 

Results of models illustrating the effect of rj and mass- 
loading on the expected cluster temperature are shown in 
Fig. and show that, not surprisingly, both mass-loading 
and rj < 1 lead to lower cluster temperatures, but also, often, 
higher X-ray luminosities. 

2.3 X-ray Luminosity Scaling Relations 

Given that the gas temperature within the cluster is compar- 
atively constant (see Fig.0 and hot (for V* > 1000 km s _1 
the temperatures are > 10 7 K, so that the bremsstrahlung 
region of the cooling curve is appropriate, i.e. A(T) oc T 1 / 2 ) 
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Figure 2. The variation of the expected X-ray temperature of 
the cluster core versus the wind velocity of the stellar component. 
Shown are model results for the cluster central temperature and 
cluster luminosity for varying r\ and the contribution of cold gas. 
The standard cluster model has M t = 10 — 4 Mq yr — 1 and R c = 
lpc. The models represented by open triangles and marked with a 
solid line are for r\ = 1 and V* = 500, 1000, 2000 and 3000 km s _1 
The models represented by filled triangles and marked with a 
dashed line are the same models except with r\ = 0.1. The models 
represented by open squares and marked with a dotted line are 
the same models except with M co id = M* 



Figure 3. The theoretical cluster X-ray luminosity versus the 
cluster wind scaling parameter X c i = M' 2 /(R C V*). Results arc 
shown for several different sets of models. The cluster radius is 
always lpc and V* = 2000 km s _1 . The open triangles (and solid 
line) are for the standard models, which have no mass injection 
{M co id = 0) and rj = 1, and stellar mass-loss rates of M* = 
10" 6 , 10" 5 , 10~ 4 and 10~ 3 Mq yr" 1 . The filled triangles (dashed 
line) have the same parameters except r] = 0.1. The open squares 
(dashed line) have r] = 1 and M co id = M* (and have very similar 
results to the previous models). The filled squares (dot-dashed 
line) have r] = 1 and M C old = 10M*. The crosses (dot-dot-dashed 
lines) have rj = 1 and M C old = 2 X 10 _ 5 Mq yr _1 . 



the following simple scaling relationships for the X-ray lu- 
minosity of clusters can be derived (ignoring for the time 
being the effects of r\ < 1 and mass-loading). 

The gas density within the cluster, n oc M*/(RcV*), so 
that if we define a cluster wind scaling parameter X c i such 
that 



Mj 
RcV* 



(10) 



then from eqn. © Lx oc X c i. We shall calculate X c i using 
the natural units, that is M in Mq yr _1 , R c in pc and V* 
in km s _1 , so that for a model with M* = 10 -4 Mq yr _1 , 
V, = 2000km s' 1 and R c = lpc, then I cI =5x 10~ 12 . 

At smaller velocities the deviations from the T 1 ^ 2 scal- 
ing for the cooling curve will alter things. For tempera- 
tures in the range of 10 6 ' 2 < T < 10 6 ' 5 K, then roughly 
A(T) oc T~ 6 and the corresponding expression for the vari- 
ation of the X-ray luminosity is 



L x oc 



(11) 



The scaling effect of variation in 77 on the cluster proper- 
ties can also be derived. A reduction in r\ effectively reduces 
the amount of thermal energy available to the cluster as a 
whole. For the same amount of mass-loss this will reduce the 
temperature, with T c i oc rj, but also reduce the velocity of 
the cluster wind and increase the density. 

For the bremsstrahlung dominant case (A(T) oc T 1//2 ), 
the density n oc M/(i? 2 t4) and V, 2 oc rj, so that here 
L x oc rj' - 5 . In the case of A(T) oc T -0 ' 6 , then L x oc t? -16 . 
Somewhat counter- intuitively then, a reduction in rj leads to 
a rise in X-ray luminosity, due to the dominant effect of the 



increasing density. Of course, for >)< 1 the emission will not 
be at X-ray energies and Lx will decrease. Results for the 
variation of the cluster X-ray luminosity versus this cluster 
scaling parameter are shown in Fig. |S| where models for a 
range of cluster wind scaling parameters X c i are shown. 

The effect of mass-loading is similar to the effect of 
lower values of rj - both lead to a drop in the average avail- 
able energy per particle. This is reflected in that models 
with 77 = 0.1 have almost identical results to a model with 
Mcoid = M t . Mass-loading leads to a decrease in tempera- 
ture and an increase in density. Model results are also shown 
in Fig. 13 Of note is that the fact that we get a very differ- 
ent relationship between Lx and X c i when the mass-loading 
injection is constant rather than a constant multiple of M* ■ 

Disentangling the combined action of lower values of rj 
and mass-loading in an individual cluster may be difficult, 
but in a statistical sample the relative influence of both 
effects may be apparent in star cluster X-ray luminosity- 
temperature plots. This is illustrated in Fig. [1] where for 
these simple models, different slopes in the cluster X-ray 
luminosity-temperature relationship are found when the ef- 
fects of rj and mass-loading are isolated and varied in a sys- 
tematic manner. Whether in reality the two effects can be 
separated so cleanly is unclear. These models have also as- 
sumed a single cluster core radius. More realistic models, 
with a range of parameters will show more scatter. 
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Table 1. The observed and predicted properties of the cluster wind from the 5 star clusters under consideration. The values for M and 
V* are obtained using eg ns . 131 and IB1 and are described more in the text. The predicted temperature, fcTrj, is calculated using eqn. l9l The 
expression for the cluster wind scaling parameter, X c i, is given by eqn. G2I The sources for observed values for Lx and kTx for each 
cluster are described in the text too. 
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Figure 4. The theoretical stellar cluster X-ray luminosity- 
temperature relationship. Shown arc model results for the mean 
cluster temperature and cluster luminosity for varying n and 
the contribution of cold gas. The standard cluster model has 
M, = 1O" 4 M yr" 1 , V* = 2000 km s" 1 and R c = lpc. The 
models represented by open triangles and marked with a solid 
line are for varying rj, with models with rj = 1.0,0.8,0.6,0.4,0.2 
and 0.1, with decreasing temperature corresponding to decreas- 
ing rj. The models represented by filled triangles and marked with 
a dashed line are for variations with the amount of cold mate- 
rial injected M co id and have rj = 1.0. The model results are for 
Mcold = 10~ 6 , 10 -5 , 10~ 4 and 10~ 3 Mq yr" 1 , with decreasing 
temperature corresponding to increasing M co id- 



3 THE SAMPLE OF NEARBY STAR 
CLUSTERS 

SSCs and other compact star clusters are relatively small 
objects with a typical size of a few pc, so even with the 
arcsecond resolution of Chandra only those that are nearby 
can be studied in any detail. In this study Chandra data 
for R136, NGC 3603, NGC 346, the Rosette Nebula and the 
Arches cluster are used. The definition of an SSC is rather 
vague and certainly the Rosette Nebula does not qualify. 
However, although it is a much smaller cluster the physics 
should be the same and so it is included. In this section we 
describe both the relevant properties of the stellar compo- 
nents of the cluster and the diffuse X-ray properties of the 



X-ray emission as seen with the Chandra satellite. In some 
cases results from the literature are used, while in others an 
analysis of the data has been performed. 

3.1 NGC 3603 

NGC 3603 is the most massive visible HII region in our galaxy, 
and is commonly regarded as a Galactic analogue for extragalactic 
SSCs. It contains a very dense concentration of stars of ages of 
~ 2 — 3Myr, including 3 WR stars in its core. In many ways it 
is similar to R136 in 30 Dor but without the surrounding cluster 
halo (Moffat et al. 1994). We assume a distance of 7 kpc. 

The characteristics of the massive stars in NGC 3603 have 
been studied by Crowther & Dessart (1998), and using this data 
and Kudritzki & Puis (2000), the mass-loss rates and terminal 
velocity data for 44 stars within 2.8pc of the centre of the cluster 
have been tabulated. These stars will dominate the stellar mass 
and energy loss into the cluster. Based on these results, the total 
stellar mass-loss rate for NGC 3603 is estimated as Af* = 2.3 X 
10 -4 Mq yr -1 and the mean weighted terminal velocity as V* = 
2844 km s _1 . The cluster core radius is taken to be 2.8pc (Moffat 
et al. 1994). 

Results from a 50ksec ACIS-I Chandra observation of 
NGC 3603 have already been presented in Moffat et al. (2002), 
and these data show strong evidence of a diffuse thermal compo- 
nent, probably associated with a cluster wind. To determine the 
spectral characteristics of the diffuse emission, point sources are 
excluded and the diffuse spectrum extracted. The spectrum has 
been fitted with an absorbed mekal model (wabs*mekal), with 
N H = 0.67 X 10 22 cm" 2 , kT x = 3.7keV, and the absorption 
corrected luminosity is Lx = 2.2 X 10 34 erg s _1 (see Table ITt . 



3.2 R136 in 30 Doradus 

R136 is the central object of 30 Doradus in the LMC, and is re- 
garded as the closest example of an SSC, and we assume a distance 
of 50 kpc (cf. Feigelson 2001). The spectral types and mass-loss 
rates of 39 stars in R136 were obtained from a study by Crowther 
& Dessart (1998), and the corresponding wind velocities from Ku- 
dritzki & Puis (2000). The cluster mass- loss rate is estimated to 
be M, = 2.6 X 10 -4 Mq yr -1 and the mean weighted terminal 
velocity as V* = 2110 km s — 1 . The cluster core radius is assumed 
to be 2pc (Crowther & Dessart 1998). 

The X-ray data of R136 was obtained from the Chandra 
archive and consists of a 28ksec ACIS-I observation of 30 Do- 
radus. The point sources within a region of radius 16" centred on 
the cluster were removed as for NGC 3603, and a spectrum was 
extracted from the remaining emission, and fitted again using an 
absorbed mekal model. The best fit spectrum to this emission 
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was an absorbed mekal model, with Njj = 0.42 X 10 22 cm 
kTx = 2.1keV, and absorption corrected Lx 
(see Table QJ. 



3.3 NGC346 

NGC 346 is the largest star formation region of the SMC, and 
contains the majority of the O stars in the SMC (Massey, Parker 
& Garmany 1989). We assume a a distance of 59kpc (Mathew- 
son, Ford & Visvanathan 1986). The spectral types and mag- 
nitudes of the stars in NGC 346 were obtained from Massey et 
al. (1989). Using a reddening correction of E(B — V) = 0.14, 
mass-loss rates were calculated from the expression in Howarth 
& Prinja (1989), accounting for the effects of lower mctallic- 
ity on stellar mass-loss rates, with M(Z) oc (Z/Zq) 0,5 (cf Ku- 
dritzki & Puis 2000 and using a value for the SMC metallicity 
of Zq/5). The cluster mass-loss rate for NGC 346 is estimated to 
be M, = 1.7 X 10" 5 M yr" 1 and the mean weighted terminal 
velocity as V, = 2282 km s _1 . The cluster core radius is taken to 
be 2pc (Massey et al. 1989). 

Chandra observed NGC 346 for a total of lOOksec in 
May 2001. The cluster lies close to a chip gaps of the ACIS-I 
instrument, but diffuse emission from the cluster is still de- 
tected. A more detailed analysis of the Chandra data of the 
NGC 346 region can be found in Naze et al. (2002a; 2002b). 
The extracted spectrum was best fit with an absorbed mekal 
model, with N H = 0.42 X 10 22 cm" 2 , kT x = l.OkeV, and 
Lx = 1-5 X 10 34 erg s _1 (absorption corrected, see Table HI . 



3.4 The Rosette Nebula, NGC 2237 

The Rosette Nebula (NGC 2237) contains an open cluster, 
NGC 2244, and is one of the more massive diffuse nebulae. It 
is situated at a distance of 1.5 kpc (Gregorio-Hetem et al. 1998). 
The cluster wind parameters of the Rosette Nebula were calcu- 
lated using the properties of the 26 O, B and A stars lying within 
4.6pc of the cluster centre, which corresponds to the rough extent 
of the X-ray emission. 

Spectral types were obtained from SIMBAD, and the corre- 
sponding wind velocities from Kudritzki & Puis (2000). Mass-loss 
rates were evaluated using the relationship of Howarth & Prinja 
(1989), assuming E(B - V) = 0.49 (Massey, Johnson, DeGioia- 
Eastwood 1995). The resulting values for the Rosette Nebula are 
M* = 2.5 X 10 -6 Mq yr -1 and V* = 2173 km s _1 . The cluster 
core radius is assumed to be 4.6pc (see above). 

Chandra observations of this region found a large number 
of point sources. On subtraction of these point sources, a weak 
diffuse component was discovered (see Montmerle et al. 2002; 
Townsley et al. 2002, in preparation). Fitting this with an 
absorbed Raymond-Smith model, the best fit parameters for the 
diffuse emission were Nh = 7 X 10 21 cm -2 , kTx = 0.6keV and 
the absorption corrected luminosity was L x = 8 X 10 32 erg s _1 
in the 0.5 — 2keV waveband. 



3.5 The Arches Cluster 

The Arches cluster is located at a projected distance of ~ 50pc 
from the Galactic centre (Cotera et al. 1996; Portegies-Zwart 
2001), and contains ~ 120 stars which have masses greater than 
20 Mq (Serabyn, Shupe & Figer 1998). It is one of the densest star 
clusters known in the local group. Blum et al. (2001) have pre- 
sented a 2/xm study of the stellar population of the Arches cluster, 
and inferred an age for the cluster of between 2 and 4.5Myr. We 



assume that the cluster core radius is 0.2pc and that the cluster 
is at a distance of 8kpc (Figer et al. 1999b). 

Due to the lack of comprehensive data on the stars in the 
Arches cluster, we have adopted the simpler approach of calcu- 
lating M* and V* using the Starburst99 models 1 . Based on the 
results of Figer et al. (1999a), for the Arches cluster we assume a 
top-heavy IMF (with a = 1.6), solar abundances, high mass-loss 
stellar evolution tracks and a total stellar mass of 1.1 X 1O 4 M 
(for stars in the mass range 1 — IOOMq), and the appropriate 
values are taken for values in the range t ~ 2 — 4.5Myr. For 
times >3.5Myr mass and energy injection from SN begin to play 
a role. For the Arches cluster, in this time span the derived 
values of the mean weighted terminal velocity are in the range 
V* = 1800 — 2810 km s — 1 , and the total stellar mass-loss is in the 
range M„ = (2.0 — 9.5) X 10 — 4 Mq yr — 1 . The corresponding large 
range in the cluster wind parameter is X c i = (0.85 — 25.0) X 10 — 10 . 

From the Starburst99 simulations it is interesting to note 
that the value of V - * > 1500 km s _1 for all times of relevance 
to young clusters. In the Starburst99 models only single stars 
are included and the energy injection rate from SN tails off at 
t ~ 35Myr (in this model). In models including binaries, mass 
and energy injection from SN will continue to longer times due 
to the processes of mass transfer and binary evolution. 

For convenience, we will adopt the values at a time of 4Myr 
for the Arches cluster, with V» = 2810 km s — , and M* = 
7.3 X 10- 4 M Q yr" 1 . These values are of course rather uncer- 
tain, and but for instance we note that direct radio detections 
of a number of the stars in the Arches cluster by Lang, Goss & 
Rodriguez (2001) might suggest a broadly similar value for the 
total mass-loss rate from the cluster. However, in their simulation 
of the X-ray properties of the Arches cluster, Raga et al. (2001) 
assumed a somewhat larger integrated mass-loss rate (60 stars 
each with an individual mass-loss rate of 10 _4 Mq yr -1 ) but a 
lower mean terminal velocity (the wind of each star was assumed 
to have a terminal velocity of 1500km s~ x ). The calculations of 
Raga et al. (2001) did not include additional mass-loading and 
in the following section we shall see that to reproduce the X-ray 
properties of the Arches (and other) clusters we find that we do 
need to include mass- loading (or a lower n), which has a sim- 
ilar effect as an enhanced mass-loss rate and reduced terminal 
velocity. 

Results from a 51ksec ACIS-I Chandra observation of the 
Arches cluster have been presented by Yusef-Zadeh et al. (2002) 
and show three main emission regions probably associated 
with the Arches cluster (designated Al— A3 by Yusef-Zadeh 
et al. 2002). Two of these components are compact (Al and 
A2) and one coincides with the core of the cluster (Al). The 
third component (A3) is more extended and part of it may 
be due to the cluster wind and part may be associated with 
fluorescent emission from a molecular cloud (C. Lang, private 
communication). Yusef-Zadeh et al. (2002) fitted the spectrum 
of this component with an additional line feature at 6.4keV. The 
total X-ray luminosity (0.2 — lOkeV) from the 3 components 
is 5 X 10 35 erg s — 1 , with component Al being the most lumi- 
nous. The spectra are fitted with a 2 component model, with 
kT\ ~ 0.8keV and kT2 ~ 6.4keV for the Al component, and we 
shall use both these temperature values in the analysis, but only 
the overall luminosity. The two compact X-ray sources could be 
genuine point sources, associated with colliding stellar winds, 
or at least one of them could also be the unresolved core of a 
cluster wind. 
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Figure 5. Left: The fitted X-ray temperature of the diffuse X-ray emission of the sample of star clusters versus the predicted X-ray 
temperature calculated from the model described in § 2. The solid line depicts equality, that is extremely efficient conversion of the 
kinetic energy in the winds of the component stars to thermal energy. Two values of observed temperature for the Arches cluster are 
plotted (see text for details). Right: The observed X-ray luminosity Lx versus the theoretical cluster wind parameter X c i feqn llQI for 
the sample clusters. Also plotted arc two results from two models shown in Fig. |3] namely models with r\ = 1 and M co id = M* (dotted 
line) and M co id = 10M* (dot-dashed line). 



4 RESULTS: COMPARISON OF THEORY AND 
OBSERVATION 

From the theory developed in § 2 and the observational results 
that we have for the small number of clusters, presented in § 3 we 
can now make some comparisons. In Table [l] we summarise the 
theoretical and observational data for the sample clusters. 

Fig. 13 (left panel) shows a plot of the observed X-ray tem- 
peratures from the Chandra observations versus the theoretically 
predicted wind temperatures for each cluster (specifically, values 
of To are plotted, calculated using can. It is clear that in all 
cases bar one the observed temperature is significantly lower than 
that the value predicted from the stellar population in the cluster. 
The exception is the case of the hotter component of the Arches 
cluster, where the discrepancy is much less. In the case of the 
Arches though, it should be remembered that the best-fit model 
to the X-ray spectrum had two components and the 2nd com- 
ponent is also plotted and is much more in line with the other 
clusters. In fact, the cooler component (kT\ ~ O.SkeV) is more 
luminous than the hotter component and it may have been bet- 
ter to use that value. Also, R136 is rather more distant than the 
rest and the possibility of unresolved point sources contaminat- 
ing the spectrum is more. If we were to ignore both the Arches 
cluster and R136 (of course, leaving only 3 data points) then a 
more coherent picture may be emerging of consistently lower ob- 
served temperatures. As has been noted before, such a reduction 
in observed temperature is a natural consequence of r) < 1 or that 
significant mass-loading is occurring. 

In § 2 the cluster wind model was also used to make pre- 
dictions concerning the scaling of the cluster X-ray luminosity 
versus the cluster wind scaling parameter X c i (see eqn. 1101 . The 
observed correlation between the measured Lx and X c i is shown 
in the right panel of Fig. [5] The most notable feature is that 
compared to the standard model in Fig. [3] the clusters are all 
overluminous by over an order of magnitude. The second feature 
is that there does seem to be a clear correlation between Lx and 
Xd, not necessarily linear, but apparently monotonic. For com- 
parison, results from model calculations, already shown in Fig. 151 
illustrate that the effects of mass-loading can roughly reproduce 
the X-ray luminosities of the cluster, with M co id = 1 — lOAif*. 



Similarly, very low values of n could yield the same result (with 
r\ ~ 0.01), though this would be in conflict with the observed 
temperatures. 

It could be argued that this correlation is partly a conse- 
quence of the simple fact that we would expect more massive 
clusters to be more X-ray luminous, and undoubtedly the value 
of X c i is most dependent on the total stellar mass-loss rate of 
the cluster (M„ ) and this will in turn scale with stellar mass. 
Larger clusters may also have larger populations of unresolved 
point sources making them appear more X-ray luminous. With 
the very limited sample we have it is difficult to go further than 
merely noting that there is a correlation between Lx and X c i. 
The Rosette nebula appears to lie at a somewhat higher luminos- 
ity than its value of X c i would imply. This could imply a larger 
comparative fraction of mass-loading as compared to the other 
clusters (see Fig.[3J. The low observed temperature, as compared 
to the expected value, would tally with this. 

In Fig. |21 we postulated that it may be possible to dis- 
criminate between the effects of a low thermalization efficiency 
and mass-loading by means of a stellar cluster X-ray luminosity 
temperature diagram. In Fig. [S] we show a plot of the observed 
values for the cluster sample. The model predicted that the 
luminosity and cluster temperature should be inversely correlated 
(though with different slopes depending on whether mass-loading 
or low thermalization efficiency was dominant). From Fig. 151 wc 
can see clearly that the data suggest that X-ray temperature and 
luminosity are probably not well correlated. 



5 SUMMARY AND CONCLUSIONS 

In this paper we have presented a comparison between the theory 
of outflows from young star clusters in a cluster wind and X-ray 
observations of diffuse X-ray emission from these clusters. The 
high spatial resolution of Chandra is necessary to begin to disen- 
tangle the diffuse emission from point source emission. However, 
even with Chandra problems remain as to being sure that the 
diffuse emission is genuinely diffuse and associated with the clus- 
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Figure 6. A stellar cluster X-ray luminosity-temperature dia- 
gram for the observed clusters, see text for more details. Note 
that the two values plotted for the Arches cluster correspond to 
the two fitted components (see text for details). 



ter wind. The major problem remains the extremely small sample 
size. More observations will alleviate this problem to some extent, 
but the objects discussed here represent the best examples of such 
a cluster wind. 

This analysis has thrown up some interesting results, the dif- 
fuse X-ray luminosity of the star clusters in the sample is corre- 
lated with the cluster wind scaling parameter X c i = / (R c V, ) 
as predicted, but that the observed X-ray temperature is not well 
correlated with the predicted X-ray temperature (though a better 
correlation may be masked by the nature of the data). It is also 
clear that from these data it is very unclear as to what is going 
on in the clusters; is there a low thermalization efficiency or is 
mass-loading the dominant mechanism? 

In this paper we have put together the results from several 
Chandra observations of nearby massive young star clusters to 
investigate the cluster wind. It is clear that many more observa- 
tions will be required to take the study further, to make a clearer 
connection between the theory and observations of cluster winds. 
Such observations will be necessary to really understand the im- 
pact that star clusters will have on the interstellar medium of 
galaxies of all shapes and sizes. 
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